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Abstract 

Experimental data from an experiment on drift-waves in plasma is presented. The 
experiment provides a space-time diagnostic and has a control parameter that per- 
mits the study of the transition from a stable plasma to a turbulent plasma. The 
biorthogonal decomposition is used to analyse the data. We introduce the notion of 
complex modulation for two-dimensional systems. We decompose the real physical 
system into complex modulated monochromatic travelling waves and give a simple 
model describing the speed doubling observed in the data as the control parameter 
increases. 
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1 introduction 



The stability of spatially extended dynamical systems and the transition to- 
wards statio-temporal chaos is still a wide open problem. Considerable progress 
has been made recently by models for pattern formation [1] and advanced 
three-wave interaction schemes [2,3]. The plasma turbulence phenomenon [4] 
is considered to be one possible application to the concept of spatio-temporal 
chaos. In particular, drift wave turbulence is of special interest for the under- 
standing of anomalous transport in magnetically confined plasmas [5]. For the 
investigation of the transition from a stable plasma to drift wave turbulence, 
the study of nonlinear model descriptions is unambiguous [6]. 
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In the present paper, we present a simple dynamical model for spatio-temporal 
bifurcations towards chaos. It is obtained from the numerical analysis of new 
experimental observations of the space-time structure of current- driven drift 
waves in cylindrical geometry. 

In Section 2, drift waves are briefly introduced. The experimental arrangement 
as well as the diagnostic tools for the measurement of the spatio-temporal 
structure of regular and turbulent drift waves are described. For the present 
investigation, three paradigmatic data sets are considered for three different 
values of the control parameter : one single propagating mode, nonlinear in- 
teraction of drift modes, and fully developed turbulence. These different dy- 
namical states result after successive bifurcations caused by an increase of 
an appropriately chosen control parameter. The powerful tool of biorthogo- 
nal decomposition [7,8] is used to analyse the experimental data. The results 
of this analysis provide the basis for the model description. In Section 3, the 
complex modulations of two-dimensional spatio-temporal systems are defined 
and the importance of modulated monochromatic travelling waves (modulated 
MTW) is underlined. The modulated MTWs are identified in Section 4 and 
characterized in Section 5. A speed doubling phenomenon is observed. A low- 
dimensional model is provided in Section 6. 



2 Experimental investigations 

Drift waves are universal instabilities of magnetically confined plasmas. In 
the present study we consider a cylindrical plasma column with a constant 
axial homogenous magnetic field B = Bz and a radial electron density profile 
ne(r). A fiuid description [9] of the plasma shows that the electrons drift in 
the azimuthal direction 9 with diamagnetic drift velocity 

v,i.^-'^L-'e (1) 

where L^^ = describes the inverse gradient length of the density 

profile. From the thermodynamic point of view, the diamagnetic drift pro- 
vides a source of free energy for the drift instability [10]. The drift instability 
propagates as a plasma wave in the azimuthal direction with diamagnetic ve- 
locity. The wave number is predominantly in the ^-direction, but has a small 
^-component to allow the electrons to flow freely along the magnetic field lines. 
The frequency of the instability is roughly given by uj* = {'m/rQ)v^ig,, where 
Tq is the position of the maximum density gradient and m is the azimuthal 
mode number. The linear analysis of a cylindrical, weakly ionized plasma has 
revealed that the growth rate of drift instabilities is strongly enhanced by 
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electron neutral collisions as well as by an additional electron drift along the 
z-axis [11]. Generally, only the onset regime of growing drift instability allows 
a linear description. In the nonlinear regime, a large variety of new dynamical 
phenomena become important (for a detailed review see Refs. [5,6]). Of par- 
ticular interest is the transition from a stable plasma to turbulence, because 
fluctuation-induced transport plays a decisive role in the plasma edge physics 
of fusion devices [12,13]. 

The experimental investigations are carried out in the central section of a large 
magnetized triple plasma device. A schematic drawing is shown in Fig. 1. A 

Fig. 1. Schematic representation of the large triple plasma device. The plasma is 
produced in two independent source chambers. Drift waves are observed in the 
plasma column in the magnetized midsection. 

homogenous quiet plasma is produced in two independent multidipole-confined 
discharge chambers [14]. For the present studies, argon is used as the filling 
gas and the degree of ionization is below 0.1%. The magnetized central section 
is separated from the plasma chambers by two electrically isolated mesh grids 
(transparency 50%). The plasma produced in the discharge chambers enters 
the central section by (i) diffusion and (ii) drift. The diffusion process (i) 
is based on the axial density gradient between the plasma chamber and the 
central section. The effect of magnetic mapping by the fringing magnetic field 
lines at the end of the central section is minimized by a particular magnetic 
compensation technique [15]. An additional electron drift (ii) is superimposed 
by biasing the grid positively with respect to the plasma potential in the source 
chamber. This electron drift is known to destabilize drift waves (see above) 
and is consequently an appropriate control parameter for the dynamics of the 
system. Experimentally, the axial electron drift is varied by the bias of the 
injection grid Uqi (cf. Fig. 1). The second grid is considered as a plasma loss 
surface and remains at anode potential, i.e : Ug2 = 0. In order to be close 
to the threshold value for the onset of the drift instability, only one plasma 
chamber is operated and a gradient-driven electron drift is present. The most 
important discharge parameters and the plasma parameters for the present 
measurements are summarized in Table 1. The radial profiles of the electron- 
density, plasma potential, and electron-temperature are plotted in Fig. 2. A 
more detailed description of the experiment can be found in Ref. [16]. The 

Fig. 2. Radial profiles of a) the electron density, b) the plasma potential, c) the 
electron temperature 

spatiotemporal structure of regular and turbulent drift waves is investigated by 
an azimuthally arranged multi-channel Langmuir probe array [16]. Each single 
Langmuir probe provides the temporal fluctuations of the plasma density. 
For the present experiment, a circular arrangement of iV = 64 probes at 
constant radial and axial position is used. The fixed radial position of the 
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Table 1 

List of discharge parameters and plasma parameters for which the experimental 
investagation is performed. 



parameter 


symbol 


value 


magnetic field 


B 


70mT 


neutral gas pressure (argon) 


v 


5.6 • 10"2Pa 


discharge voltage 


Ua 


60V 


discharge current 


Id 


13A 


plasma radius (FWHM) 


R 


0.1m 


plasma length 


L 


1.8m 


e drift velocity 


Vd,c 


< OAvth,e 


electron density (center) 


ric 


5 • lO^Ocm-3 


electron temperature (center) 




1.5eV 


ion temperature (center) 


Ti 


0.03eV 



probes, imposed by technical constraints, doesn't allow to address here the 
interesting problem of the radial profile of the perturbation. The probe array 
provides the temporal evolution of the spatial structure of drift waves which 
propagate mainly in the azimuthal direction. The temporal resolution is given 
by the maximum sample rate of the aquisition system {At = l/xs), and the 
spatial resolution is given by the azimuthal angle between each two probes 
(Ax = 27r/64). The data is stored as an x M-matrix Uij = n^(iAx^jAt) 
where = 64 (space) and M = 2048 (time). Note that investigations of 
spatiotemporal phenomena can be done with one probe by a method based 
on conditional averaging [17]. 

In order to study the bifurcation behaviour of drift waves, the data sets of three 
representative dynamical states are considered: Ui a single monochromatic 
drift mode, U2 drift modes with nonlinear interaction, U3 turbulent drift waves. 
The different dynamical states Uk are recorded successively by increasing the 
accessable control parameter C/qi- In Fig- 3, the three data sets Uk are shown 
as greylevel plots. Data set Ui [Fig. 3(a)] is a single propagating mode, data 

Fig. 3. Greylevel plots of the three drift wave data sets. Grey corresponds to the 
zero level, black is the minimum, and white is the maximum deviation from the 
equlibrium electron density in the plasma. The vertical axis is labelled with the 
probe number index, the horizontal axis is the time in multiples of At = l^us. 
Control parameter values arc (a) Ugi = 4.5V, (b) Ugi = 4.6V, (c) Ugi = lOV. 

set U2 [Fig. 3(b)] shows interacting drift modes, and data set W3 [Fig. 3(c)] is 
the state of strong drift wave turbulence. These three data sets are the basis 
for the following analysis and model description. For this purpose the plasma 
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is considered as a spatiotemporal dynamical system, whose state is described 
by a function u^{x,t) where e represents the complete set of experimental 
parameters. For the present investigations only one parameter, the grid bias 
Uqi, is varied. The grid bias has been chosen as control parameter because 
it allows to consider that all other parameters remain in good approximation 
constant. 

Each dynamical state given by the data set Uk corresponds to a control pa- 
rameter value ejfc. 

The analysis tool used here is the biorthogonal decomposition (BOD). This 
tool provides a way to study in the space and time properties simultaneously. 
We present here just the most important features of the BOD. (For more 
details see Ref. [7]) 

Suppose that our system is described by a function u{x, t) defined on a spatial 

range X and a temporal interval T. In the experimental situation, X is the 
domain of the azimuthal angle x in cylindrical coordinates, i.e. X = [0,27r]. 
The biorthogonal decomposition provides the smallest linear subspace x(X) 
containing the phase space trajectory (described as time t runs) defined by 



The set of all the vectors is the trajectory and the evolution of is the 
dynamics of the system. 

The biorthogonal decomposition also provides the smallest linear subspace 
x(T') containing the spatial structure (described as the spatial position x 
varies) defined by : 



In the present paper, the scalar product is used to define H[X) and H(T), 
the Hilbert spaces of the functions of x defined on X, and the functions of t 
defined on T respectively. The BOD is the spectral analysis of the operator 
[/, which acts from H{X) into H{T) : 



where U defines a one-to-one relation between the vectors of x(X) and x(T), 
the orthogonal complements of the kernels of U and its adjoint U*. 



\/x e X, ^t{^) = u{x,t). 



(2) 



VteT, ^^{t) = u{x,t). 



(3) 




(4) 
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If U is compact, as it is here the BOD decomposes u{x,t) into temporal and 
spatial orthogonal modes and u{x,t) can be written as follows : 

u{x,t) = ^a„^/'„(t)0„(a;), (5) 

with ai > a2 > . . . > 0, and the orthogonality relations (pm) = ^n,m and 
i'ipn, i^m) — Sn,m- The (pn are called topos, and the ipn chronos. 



3 Modulated travelling waves 

In this section we discuss the conditions for a two-dimensional system to 
be considered a modulated travelling wave. The introduction of modulated 
travelling waves is of interest as we shall see in the drift wave turbulence, 
where the decomposition of a high-dimensional system into a set of two- 
dimensional sub-systems allows us to focus on the dominating properties of 
the dynamical behaviour. The way in which a monochromatic travelling wave 
may be deformed in our case of study is now analysed in the following theorm 
for which we need first to introduce some simple definitions concerning a 
complex formalism for the corresponding modulations. 

Let us now specify what is meant by two-dimensional projections, of a N- 
dimensional system be described by its BOD 

N 

u{x, t) = Y^ akipk{t)(t)k{x)- (6) 

k=l 

defined on the space interval X — [0, 27r] and on the time interval T — [to, ii]. 

Definition 1 The projection of a system u{x,t) onto two vectors of index 
m and n is the two dimensional system Um,n{x,t) described by : 

t) = am1prnit)(pmix) + an1pn{t)(pn{x) (7) 

Thus the projection of the dynamics (associated with the spatio-temporal 
system u{x, t)) onto the eigenvectors of index m and n is simply the projection 
of onto the two topos 0^ and 0„, i.e. 

yxeX, er»=tim,n(2:,t). (8) 

The projection of the spatial structure (also given by the spatio-temporal 
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Fig. 4. The operator U maps the subspace of xi^) spanned by (pi and 02 into the 
subspace of x(T') spanned by ■^i and '02- 

system u{x, t)) onto the eigenvectors of index m and n is simply the projection 
of onto the two chronos ■0m and ipn- 

VteT, C"(t)=«„,„(a:,t) (9) 

The schematic drawing Fig.4 illustrates this projection process of the dynamics 
and the spatial structure. 

Definition 2 A modulatrix A4 is a pair of complex valued continuous func- 
tions M{x) and N{t). M is called the spatial modulatrix, and N the temporal 
modulatrix. 

Note that each complex-valued continuous function is the parametric repre- 
sentation of an arc. Using the vocabulary of complex analysis (see for instance 
[18]), the arc Fx is defined as 

Fx : M{x), 0<x<2n. (10) 

Definition 3 A modulatrix M. = (M, N) is a continuous phase modula- 
trix if M and N can be written as 

M(x) = A(x)e'-^^^) (11) 

and 

N{t) = B{t)e'^^'\ (12) 

where A, B, F , and G are continuous functions. It will he called a regular 
continuous phase modulatrix if moreover the increase of the argument 
of F , Dp = F{2t{) — F{fS), and the increase of the argument of G, Dq — 
G(ti) — G(to), are both equal to zero. 

It is known that a continuous complex-valued function can be written as a 
product of a continuous modulus function and a continuous argument function 
if the function is never equal to (see [18]). The modulus function is unique 
and two different argument functions differ by a constant integral multiple 
of 2n. The real numbers Dp = F{2n) - F(0) and Dq = G{ti) - G{to) are 
independent of the choice of the argument functions F and G [18]. 

Definition 4 The spatial complexification of two-dimensional system u{x, t) 
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whose BOD is 

u{x, t) = aiilJi{t)(j)i{x) + a2ilJ2{t)(l)2{x) 



(13) 



is 

Z{x) = ai4)i{x) + ia2<i)2{x)- (14) 

The corresponding temporal complexification is 

y(t) =ai^i(t) +m2i^'2(t)- (15) 

Note that the spatial complexification is a representation of the spatial struc- 
ture in the complex plane, and the temporal complexification is a representa- 
tion of the dynamics in the complex plane. 

If the complexifications are never zero, the argument function can be intro- 
duced by the following definitions : 

Definition 5 A spatial complexification Z{x) is a phase continuous com- 
plexification if it can be written as : 

Z(x) = C{x)e"^^^\ (16) 

Definition 6 A temporal complexification Y (t) is a phase continuous com- 
plexification if it can be written as : 

Y{t) = D{t)e''^^'\ (17) 

Of particular interest is the case in which X is the circle. 

Definition 7 Let u{x,t) be a two-dimensional system defined on a circle X. 

Its spatial complexification Z{x) is assumed to be always non-zero and contin- 
uous on X . Then Z{x) can be written as Z{x) = C{x)e'^'^^^\ As Z{0) = Z(2'k) 
we have Q{2tt) = Q{0) + n27r where n is an integer. The number n is called 
the spatial winding number of the system. 

Let us now introduce the notion of modulation of a system : 

Definition 8 Let u{x,t) be a two-dimensional system which BOD is given 
by 13. Let M. — {M{x), N(t)) be a (regular) modulatrix. A spatiotemporal 
(regular) modulation ofu{x,t) is the system u'{x,t) defined by 

u'{x,t) = a'Mt)<j>[{x) + a'2^lJ'2{t)cj>'2{x) (18) 
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with 



a[(f)[{x) + iQ;202(^) = M{x){ai(f)i{x) + ic(2(j)2{x)) (19) 

and 

a[ip[{t) + ia'^ip'^it) = N{t){aiipi{t) + ia2Mt))- (20) 

Note that in general 4>'i, ip2, and 02 are not eigenvectors of u'. 

Let us investigate the effect of a modulation on the dimension of the system. 

Theorem 9 Let u{x,t) be a system described by two phase continuous com- 
plexifications Z{x) = C(a;)e^^(^) and Y{t) = D{t)e'^^^\ Let M = {M,N) be a 
continuous phase modulatrix with M{x) = y4(a;)e*^^^^ and N{t) — B{t)e^^^*\ 
The dimension of the modulated system is reduced to one if 

Vx e X, Q{x) + F{x) - zi (21) 

or 

Vt e T, R{t) + G{t) = Z2 (22) 

where Zi and Z2 are two complex numbers. Otherwise the dimension of the 
modulated system remains equal to two. 

PROOF. Let us first apply the spatial modulation. Let Z'{x) be the spatial 
complexification of the spatialy modulated system u'{x,t) : 

Z'{x) = M{x)Z{x) (23) 

The dimension of the spatial structure is reduced only if it is embedded in a 
segment, i.e., if Z' can be written as Z'{x) — C"(x)e'*i where qi is real. In this 
case, 

u'{x,t) = C (x) cos qiijjiit) + C'{x) sin qi'ilj2{t), (24) 

it is clear that the new system is one-dimensional, more precisely 

u'{x,t) = a[ilj[{t)(p[{x) (25) 
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with q;^ = 01 = C /a'-^, and = cosgi'0i + singi'02- 

However, 

C\x)e"'' = yl(x)e'-^(^)C(x)e'^(^) (26) 

This equation imphes G X, + F[x) = 6*"^% where Zi = e*^^ Con- 

versely, if we assume that M{x) = A{x)e'^^'^'^~^^^^\ then the modulated system 
becomes one-dimensional. We can show in the same way that the reduction 
of the dimension of the dynamics after a temporal modulation is equivalent 
to (22). □ 

Theorem 10 A regular spatial modulation does not change the spatial wind- 
ing number. 

PROOF. Let Z{x) = C{x)e^^^^^ be a complexification with a winding num- 
ber n. We then have Q{2ti) = Q{Q) -\-n2Ti. Let M{x) = A{x)e^^^^^ be a regular 
spatial modulation. Therefore, because X is here the circle, F{2Ti) = F{0). 
The spatial complexification of the modulated system is by definition 

Z'{x) = C"(x)e^^'W (27) 

with C'{x) = C{x)A{x) and Q'{x) = Q{x) + F{x). Then using the values of 
F{x) and Q{x) at x = and x = 27r we get (5'(27r) = g'(0) + n27r. □ 

Note that this theorem implies that in order to change the winding number 
by modulation, a non-regular modulation is required. 

Theorem 11 Let u{x,t) be a two-dimensional system defined on the circle 
X whose spatial complexification Z{x) and its temporal complexification Y{t) 
never vanish. 

Then there exists one unique regular modulatrix and a pair consisting of an 
integer k and a real lo such that u{x, t) is the modulation of a monochromatic 
travelling wave Uo{x,t) — cos (kx) cos (cut) + sin{kx)sin{ujt) . 

PROOF. Let us write the spatial complexification Z{x) as 

Z{x) = C(x)e^'2(^) (28) 
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and the temporal complexification as 



Y{t) 




(29) 



We set 



k 



1 



(30) 




and 



1 



(31) 




where Dq = Q{27t) — Q{0) and Dr = R{ti) — R{to) are the increases of the 
arguments of Z and Y respectively. The numbers k and cu are unique, because 
Dq and Dr are unique. Because X is a circle, the winding number k is an 
integer. 

Let us define : 



The modulatrix A4 = {M,N) is regular and u{x,t) is the modulation of a 
system -^0(3;,^) by Ai. The spatial complexification and the temporal com- 
plexification of Uo{x,t) are respectively Zo{x) = e**^^ and Yo{t) = e*'^*. The 
system uo{x, t) is thus of the form 



The functions Q{x) and R(t) are defined modulo 2%, so C{x) and D{t) are 
unique. Thus the functions M and N are unique. □ 

This theorem provides a way to consider a two-dimensional experimental sys- 
tem as a spatio-temporal modulation of a monochromatic travelling wave. In 
particular, a continuous deformation of a wave is described by a modulatrix 
M.^ = (Mg, iVg) that depends on the control parameter e. In the case where X 
is the circle, we can define the winding number, and the jump in the winding 



M{x) = C(x)e^[^(^)-'=^1 



(32) 



and 



N{t) = L>(i)e^[^(*)-^*] 



(33) 



uo{x,t) — cos{kx)cos{ijjt) + sin{kx)sin{ujt). 



(34) 
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Fig. 5. Plot of the weights for the data in y-log scale (left) and plot of the ratios 
an+i/fin with respect to n (right), (a) data I4i (b) data U2{c) data U3. They ex- 
hibit the degeneracy of the weights which may be associated to a spatio-temporal 

symmetry 

number that is observed for a certain value of the control parameter will cor- 
respond to a bifurcation. Note the analogy of the complexifications with the 
loops in the order parameter space in the study of the topological defects [19]. 

Our complex modulation is a generalisation of the real modulation introduced 
in [20] . The modulations introduced in [20] are not appropriate for the descrip- 
tion of the deformations of a function 0i(x) (resp. ipi{t)) that has a zero which 
shifts as e varies. Instead, we can describe the deformation of such a function 
provided that we find a complementary function 02 (a;) (resp. tp2{t))thai never 
vanishes for the value of x (resp. t) where (f)i{x) (resp. "01 (^)) has a zero. 
Note that in the absence of resonances, as in the case of real modulatrix 
A4 = (M, A^), the new eigenvectors are the modulation of the non perturbed 
system, as in Ref. [20]. 



4 Identifying the modulated monochromatic traveUing waves 

The aim of this section is to simplify the description of the dynamical be- 
haviour of the present system by the identification of the two-dimensional 
structures existing in the system for different values of the control parameter 
e. In order to identify two-dimensional subsystems whose spatial and tempo- 
ral complexification never vanishes, the idea that this system is a deformed 
monochromatic travelling wave (MTW) is used. More precisely, a MTW has 
two basic properties : (i) the degeneracy of the weights and (ii)the Fourier 
transform of the chronos and topos are delta functions. 

We show in Fig. 5 the plot of the weight on a logarithmic scale. In the analysis 
of the plot of weights, we have to look for pairs of eigenvalues that arc degen- 
erated. The degeneracy of pairs is obvious only for the first two eigenvalues of 
each data set. However, if the ratio a„/ a„+i is plotted versus n, the degeneracy 
of all weights can be quantified (see Fig. 5) 

We are thus looking for values of n such that a„/a„+i is close to unity. This 
plot provides a way to couple eigenvalues : two eigenvalues a„ and a„+i are 
considered to be coupled if the ratio a„/a„+i forms a local maximum in the 
plot ratios. Following the previous rule, we obtain the pairs listed in Table. 2 
that we call structures. 

Notice that, considering only the weight distribution is not sometimes suffi- 
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Table 2 

Pairs obtained after the analysis of the weigths 
data set 





1 


-2 


3- 


-4 




6- 


-7 


9-10 


12- 


13 


14- 


-15 


16-17 




1- 


-2 


3- 


-4 




7- 


-8 


9-10 


12- 


-13 


14- 


-15 


16-17 


U3 


1- 


-2 






6-7 


8- 


-9 


10-11 


12- 


-13 


14- 


-15 


16-17 



Fig. 6. Modulus squared Fourier transforms of the chronos. The three columns 
correspond to the three data sets. The x-axis is the frequency. The frequency unit 
is du = 2048.10~^s~^. The y axis is labelled by the eigenvector index. 

Fig. 7. Modulus squared Fourier transforms of the topos. The three columns cor- 
respond to the three data sets. The x-axis is the spatial frequency. The frequency 
unit is the spatial loop. The y axis is labelled by the eigenvector index. 

cient to decide the pairing. Such is, for instance, the case in data U3 for the a„, 
n = 3, 4, 5. A careful look to the corresponding chronos and topos, specially 
by considering their Fourier transform, will give, in this case an unambiguous 
pairing. 

The next step in the identification of modulated MTW is to consider pairs of 
functions that have the same spatial (temporal) frequencies. An appropriate 
way to detect such modulated sine and cosine functions is to apply a Fourier 
transform to the eigenf unctions. The plot of the square modulus of the Fourier 
transforms of the chronos and topos are represented in Fig. 6 and Fig. 7. 

A single-peaked spectrum is found for the first two chronos and topos of each 
data set. Pairs of eigenvectors are found by inspection : each two Fourier 
spectra with almost or exactely the same structure indicates the presence of a 
pair. Using the conditions on the Fourier transform, Table 2 can be improved. 
Table 3 give the new list of pairs. The pairs are labelled Wk, where k is the 
index of the wave associated with the pair of eigenvectors. 

Table 3 

Pairs obtained after the analysis of the Fourier spectrum of the eigenvectors. The 
modulated MTW associated with those pairs are labeled Wk, where k is the index 
of the wave. 



data set 


Wl 


W2 


W3 


W4^ 


W5 






Ws 


Ui 


1-2 


3-4 


5-6 


7-8 


9-10 


12 13 


14-15 


16-17 


U2 


1-2 


3-4 


5-6 


7-8 


9-10 


12-13 


14-15 


16-17 


Us 


1-2 


3-4 


5-6 


7-8 


9-10 


12-13 


14-15 


16-17 



Note that the Eigenvector 11 for the sets Ui and U2 and Eigenvector 7 for the 
set U3 do not belong to a pair. 
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Fig. 8. Spatial structures associated with the first eight modulated travelling waves, 
a) Data set Ui, b) Data set U2-, c) Data set W3. Each square subfigure is labelled in 
the corner by the index of the eigenvectors spanning the structure. 

Note also that the degeneracy of the weights a% and aj doesn't correspond to 
a modulated travelling wave. This is because the energy of the eigenfunction 
7 wich is a global oscillation of the plasma without any propagation is close 
to the energy of the third travelling wave. 



5 Organization of the different modulated MTWs 

In the previous section, two-dimensional subsystems which are modulated 
travelling waves were identified . In this section, these waves are studied with 
respect to both the eigenvector indexes and the control parameter. 

First the spatial and temporal frequencies k and uj are determined. Actually, 
Theorem 11 provides an explicit way to compute k and uj. However, this direct 
approach is cumbersome in practice. Therefore the winding number k is ob- 
tained in a different manner : The first eight spatial structures associated with 
the first eight modulated travelling waves are shown in Fig. 8. The winding 
number becomes easier to inspect in the plot of the spatial structures (as in 
the case of the Fourier transform) as the control parameter becomes larger. 
This shows that the spatial modulation becomes more regular as e increases. 

Fig. 8 shows that the winding number for the data set 1J\ is well defined only 
for the first wave wi (The labelling is defined in the Table 3). For the data set 
ZY2, the spatial winding number is defined for the three first waves, and for the 
data set W3 it is well defined for all eight waves studied. The winding numbers 
found are shown in the Table 4. The temporal frequencies u) associated with 
these waves are determined from the Fourier spectra (Fig. 6). A broad Fourier 
spectra corresponds to the presence of zeros in the temporal complexifications. 
The UJ is well-defined from the Fourier spectra only for the first wave of the 
data set Ui and U2, and for the first three waves for the data set W3. 

The phase speed of a modulated wave is defined as the ratio of the tem- 
poral frequency a; to the spatial frequency k of the associated unmodulated 
monochromatic wave. Considering the first waves given in Tab. 4 it is noted 
that a speed doubling occurs when the control parameter e increases from the 
value that belongs to the data set IA2 to that of the data set lA^. Indeed, the 
phase speed associated with the waves wi, W2, and in the data set W3 is 
very close to twice the speed of the first travelling wave of the data U2- In the 
next section, a simple model for this spatio-temporal bifurcation is discussed. 
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Table 4 

Spatial winding numbers k and temporal frequencies lo associated with the first 
eight waves. The wave number i is labelled wave Wi 
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Fig. 9. Graylevel plots of the two dimensional structures. Each column is associated 
with a data set. For each column, the letters a) . . . h) indicate the wave wi,. . . ,ws 
that is plotted. For each figure, the x-axis is time and the y-axis space. 

Another way to study the effect of the spatial modulation is to plot a graylevel 
plot of the two dimensional restriction Um,n{x,t) associated with the modu- 
lated wave. The graylevel plots for the first eight waves of each data set is 
presented in Fig. 9. 

In those plots a phase spatial modulation implies a distortion of the wave front, 
this distorsion depending only on x, and an amplitude spatial modulation 
corresponding to global crashes in the amplitude of the waves, those crashes 
having the same intensity for fixed positions. In the graylevel plot, the phase 
modulation is easier to observe than in the probe structure, where a phase 
modulation implies only a non- uniform spacing of the vectors ^x- 



The graylevel plots show that in the data sets Ui and W2, the first wave has a 
phase defect that localy makes the wave front more vertical. The other waves 
are strongly phase modulated with strong tearing of the wave fronts. However, 
the wave fronts of the waves in the data set U2 are more regular than in the 
data set Ui. In the data set ZY3, all the strong phase defects have disappeared 
or been reduced. 



In the gray-level plots the strong amplitude modulations can be observed as 
well. The global crashes of the amplitude of the wave correspond to an ampli- 
tude modulation close to zero. This amplitude modulation is itself chaotic. 

The data set W3 is thus a state where chaos is essentially present in time. The 
spatial structures are, on the contrary, fairly regular. The temporal chaos is 
then due to a chaotic modulation acting on a non-chaotic spatio-temporal 
structure, i.e. a travelling wave. Note that, furthermore, in this case we are 
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faced to an homogeneous turbulence [21] as it is seen in Fig. 6. Also note that 
the chaotic modulation of monochromatic waves has a close relationship with 
the three-wave interaction model of the drift wave instability [6]. However, 
it one may ask if the restriction to three waves is pertinent here. Indeed, our 
study reveals a set of at least eight waves relevant to the data set Uz- 

It is easier to study the spatial amplitude modulation in the probe structure 
plots [Fig.8]. The modulation of the first wave increases with e. In contrast, 
the next waves become more regular as e increases. The spatial structure of 
the data set coincides with a large scale of eigenvectors (m = 3, . . . , 19). 

To study these large scale phenomena, the plot of the weights is considered 
[Fig. 5] . Neglecting the pairwise degeneracy, the weights decrease exponen- 
tially with n in well defined regions. In particular, in the ditribution of the 
weights for the data set W3, the boundaries of such a region are given by the 
indices rii = 7 and n2 = 20. These boundaries correspond to a pronounced 
broadening of the Fourier spectrum in the low-frequency regime [cf. Fig. 6]. 
In this region, the spatial structures have a well defined winding number as 
shown in Fig. 8. Note that the spectrum of the chronos is bounded by a fre- 
quency near 60 in the data set Ui and W2. In the data set W3 this bound has 
been shifted to 250. Even in the eigenvectors with a broad spectrum, i.e. those 
with an index greater than 23, the frequencies higher than 250 are nonexis- 
tent. However, the spatial frequencies are not bounded and increase with the 
index of the eigenvector. 

Remark 12 Even if the radial dependency of the electric field is not directly 
accessible to the present diagnostic, we may wonder that such effects is present 
giving rise to a collective rotation of the whole plasma column and therefore 
to a shift in the (Fourier) dispersion relation for the waves. This effect seems 
present in the frequencies reported in Table 4, for the data Us. However, since 
we restrict ourselves to each two dimensional structures , this fact would have 
no effect in our analysis. It simply changes the position of x{T) in the space of 
the functions of time H{X) by a global rotation. We thank one of the referees 
to have pointed us this question. 



6 Model for the speed doubling 

In the previous section, it was noted that a speed doubling occurs when the 
control parameter increases. The data set IA2 corresponded to a state before 
this bifurcation and the data set corresponded to a state after the bifur- 
cation. In this section we give a simple model for this speed doubling. Our 
model is built with two two-dimensional structures. 
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Fig. 10. A spatial modulation can imply a topological change in the spatial structure, 
which corresponds to a change of the speed of the modulated travelling wave. The 
spatial structures, the trajectories, and a sketch of the crest of the wave are plotted 
for four values of the control parameter. 

Fig. 11. Model for the speed doubling. The gray level plot of the modulated MTW 
is shown for four values of e. The corresponding spatial structures are shown in the 
right column. 

We first describe one of the two structures. The figure Fig. 10 shows how a 
spatial modulation can change the winding number and the speed of the speed 
of the eigenstructure. 

The phase speed (cf. Section 5) is defined as a ratio. In order to define the 
states to model uniquely, the speed alone is not sufficient. The model structure 
describes an evolution from a state described by the function u^i{x,t) ~ 
cos(2a;) cos(cjt) + sin(2x) sin(cjt) to a state described by the function ui{x, t) — 
cos(a;) cos(ci;t) + sin(x) sin(ci;t). 

Let us consider a temporal complexification which is independant on e defined 

by 

Yo{t) = e'^K (35) 

and the corresponding spatial complexification by 

Z,{x) = (1 - e)e*2" + (1 + e)e'^ (36) 

These spatial and temporal complexifications correspond for e = — 1 to the 
complcxifications of U-i{x,t) and for e = 1 to the complexification of Ui{x,t). 
Fig. 11 shows the plots for four values of the control parameter e, the spatial 
complexification Z^{x) in the right column, and the gray level representation 
of the two dimensionnal system u^{x, t). 

The spatial complexification is never equal to zero except for the value e = 0. 
For this value, the spatial complexification vanishes for a; = tt. For negative 
e the winding number is 2, and for positive e it is 1. However, for e = 0, 
the winding number is not defined. A bifurcation in the winding number thus 
occurs. 

A regular spatial modulation is just defined (i) before the bifurcation by divid- 
ing the spatial complexification Z^{x) by Z_i{x) (ii) after the bifurcation by 
dividing the spatial complexification Z^{x) by Zi{x). The spatial modulation 
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before the bifurcation M^{x) is defined by 

M^^{x) = (1 - e) + (1 + e)e-'^ 



(37) 



M^{x) is a regular modulation because (1 — e) > (1 + e) when e < 0. The 
spatial modulation after the bifurcation M^{x) is defined by 

M^{x) = (1 + e) + (1 - e)e^^ (38) 

M^{x) is a regular modulation after the bifurcation because (1 + e) > (1 — e) 
when e < 0. The spatial defect which occurs at x = tt for e = permits the 
wave front to change its shape as shown in Fig. 11. The two (unnormahzed) 
topos are the real part and the imaginary part of the spatial complexification 

(j)^(x) = (1 - e) cos(2a;) + (1 + e) cos(a;), (39) 

02(x) = (1 - e) sin(2x) + (1 + e) sin(x). (40) 

The chronos are left unchanged 

i/jiit) ^cos{ujt), (41) 
^2(2;) = sin(a;t). (42) 

Both the spatial behaviour (see Fig 8, Fig 9) — i.e. the bifurcation of a winding 
number from the value 2 to the value 1 — and the temporal behaviour (see 
Fig 6, Fig 9) — i.e. a frequency lu roughly independant of the control parameter 
e shows that this structure corresponds to the modulated MTW Wi for the 
data sets Ui and U2 and W2 for the data sets U3. (The notations Wk were 
introduced in the Table 3). 

In the same way we model the structure corresponding to the second pair 
in Ui and U2 and the first in U^. Notice that in this case the bifurcation is 
due to the simultaneous deformation of the temporal modulation and of the 
corresponding spatial modulation. 

Finally, the model is built by glueing the two structures weighted by two cor- 
responding eigenvalues (the relative energies of the two structures). According 
to the general theory of bifurcations described by the BOD [22] , these ener- 
gies cross at the bifurcation and this is the reason why the energies of the two 
structures are interchanged when passing from U2 to U3. 
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7 Summciry and Conclusions 

The BOD of the experimental data of the drift-waves experiment showed the 
importance of the two-dimensional sub-systems. It has been shown that the 
two-dimensional structures are complex-modulated monochromatic travelling 
waves (modulated MTW). A spatial regularity counterpart of a temporal chaos 
has been discovered in the most turbulent data set. The most important fea- 
ture of the evolution of the system with the control parameter e, i.e. the speed 
doubling, has been modeled as well. The model consists of a pair of simple 
two-dimensional structures that undergo an exchange of energy as the value 
of the control parameter varies. 

A study of new experimental data will be done in the future in order to im- 
prove the model and better characterize the bifurcation. The regular behaviour 
also needs to be better understood. Furthermore, the connection between the 
modulated MTW and the waves-interactions models for the drift-waves will 
be investigate in the future. 
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